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We consider the time delay of massive, non-relativistic, one-dimensional particles due to a tun- 
neling potential. In this setting the well-known Hartman effect asserts that often the sub-ensemble 
of particles going through the tunnel seems to cross the tunnel region instantaneously. An obstacle 
to the utilization of this effect for getting faster signals is the exponential damping by the tunnel, 
so there seems to be a trade-off between speedup and intensity. In this paper we prove that this 
trade-off is never in favor of faster signals: the probability for a signal to reach its destination before 
some deadline is always reduced by the tunnel, for arbitrary incoming states, arbitrary positive and 
compactly supported tunnel potentials, and arbitrary detectors. More specifically, we show this for 
■ several different ways to define "the same incoming state" and "the same detector" when comparing 

the settings with and without tunnel potential. The arrival time measurements are expressed in the 
time-covariant approach, but we also allow the detection to be a localization measurement at a later 
£S| ■ time. 

PACS: 03.65.Db, 03.65.Nk 
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INTRODUCTION 



Questions related to the tunneling phenomenon have been actively studied since the early days of quantum me- 
chanics, and some of them are still not resolved. In the simple case of a massive particle moving in one dimension 
through a localized (tunnel) potential, the question of the "time spent in the tunnel" is especially interesting, and 
has given rise to extensive discussion (see e.g. 0, 0, [HI, Hi], Hlj and the references therein). Some difficulties in 
Q-i dealing with this problem are rooted in the absence of a selfadjoint "time operator" ("Pauli's Theorem" [22|). Instead, 
l— ~ 1 ■ one has to use more general framework of positive operator measures (POMs) @, [IB, Gil H3|- For a survey on time in 
! \ quantum mechanics, see [2C| . 

An old observation related to tunneling times is the so called Hartman effect [Tj|, which states that the transmitted 
' part of a wave function appears to move faster through the tunnel than the corresponding free state. More precisely, 
after a long rectangular barrier and for a wave function of narrow momentum distribution, in leading order the 
transmitted pulse appears at the end of the tunnel instantaneously. Therefore, it has been suggested that this effect 
means superluminal signal transport [3, El- However, all this is only true for the shape of the wave function 
t~^. ■ disregarding normalization. But obviously, especially for long tunnels, for which the gain in speed would only be 
noticeable, the transmission probability is exponentially small. So in any attempt to utilize the Hartman effect for a 
faster signal transmission, we would have to analyze the trade-off between transmission probability and transmission 
speed. 

The main result of this paper is that this trade-off is always trivial: when damping is taken into account, transmission 
through a tunnel will always slow down the signal. Figures 1 and 2 show a sketch of the result of an arrival time 




FIG. 1: Arrival probability density for free (dashed line) and transmitted (solid line) particles 
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FIG. 2: Probability for arrival up to time t for free (dashed line) and transmitted (solid line) particles 

measurement in a possible situation: The arrival time probability density for the transmitted particles peaks earlier 
than for the free particles. The density may even become larger at some times. But if we look at the integrated 
density, i.e., the probability for the particles to arrive before a given deadline t, plotted here as a function of t, then 
the free particles win. 

This is true in remarkable generality: for any incoming state, any tunnel potential, and any detector. It is even true 
for several, in general inequivalent approaches to formalizing the rules of this race. Indeed we need to choose a precise 
notion of arrival detection, of the equality of initial wave functions, and of the equality of detectors for tunneling 
and free dynamics. The equality of initial states is a non-trivial issue, because the two states are subject to different 
dynamical evolutions. So at least we need to fix a reference time. For this we have two choices, namely either a fixed 
time (set to zero by convention), or t = — oo, i.e., asymptotic equality of the incoming states in the sense of scattering 
theory. In this case we fix the direction from which the particles are coming by choosing input states with positive 
momentum. 

On the detection side, a natural choice is to describe the detectors as covariant arrival time observables dam. 
Again this raises the issue of how to compare the two cases, because the covariance condition explicitly depends on the 
time evolution, and an observable can be covariant with respect to only one of them. Again scattering theory helps, 
by defining a bijective correspondence between the respective sets of covariant observables: we identify observables, 
which give the same probability distributions on states coinciding for t — > +00. This identification is also natural 
for including finer descriptions of the detection process. For example, we could modify both time evolutions by 
including an imaginary "optical" potential, resulting in contraction semigroups rather than unitary groups. The loss 
of normalization is then interpreted as arrival probability. Even more realistically, we could model the detector by 
a system in a bound state, interacting with the particle through a potential, and getting ionized in the detection 
process. The ionization time (as measured by a covariant observable on the escaping electron) is then taken as the 
detection time. It will be shown elsewhere how these ideas lead to special cases of covariant measurements. 

In this paper we also look at another way of setting up the finish line for the race: For particles traveling in the 
positive re-direction, any time t, and any position a behind (=to the right of) the tunnel, we can replace the event "the 
particle has arrived at point a before time f" by "the particle is located in the half axis [a, 00) at time tf'. The latter 
statement requires only a position measurement, and hence does not require the theory of arrival times. Although the 
two statements are not equivalent, and correspond to different effect operators (the first probability is, by definition, 
increasing as a function oft, but the second is not in general increasing), they will be qualitatively similar, and equal 
in the classical limit. Note that both the arrival times and the localization observables (see e.g. [2j, [3j]) are represented 
by positive operator measures, constrained by a covariance condition. In order to emphasize the analogy with arrival 
times, for which a projective measurement does not exist, we also allow localization observables, which are general 
positive operator measures ("POM") 1 , rather than just the standard position observable. 

To summarize, we are looking at the following three approaches to our problem: 
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1 More commonly called "POVMs" for positive operator valued measures 
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The combination of the time-zero initialization with time-covariant measurements is conspicuously absent from this 
table, because the identification of detectors under different dynamics requires asymptotic scattering theory, which 
we wanted to avoid in this setting. 

There would be more possibilities for the detection, with more realistic detector descriptions, and some of these 
are now under investigation. Approaches I and II have been discussed in a recent paper by Werner and Ruschhaupt 
(to be published, see also a conference report by Werner at the 40th Toruh Symposium, June 2008). Here we have 
added the "time zero approach", as well as formulated the treatment in a way that brings the use of positive operator 
measures and covariance to the front. 



II. PRELIMINARIES AND NOTATIONS 

As mentioned in the introduction, we have three distinct approaches, each formulated in terms of covariant positive 
operator measures. Because of the covariance property (which we will precisely define shortly below), the mathematical 
structure of the problem will be the same in each case, once the observable is suitably transformed. More specifically, 
time observables are defined via time translations and localization observables via space translations; in both cases, 
the structure is the same in the Hilbert space where the translation generator acts multiplicatively. Accordingly, 
we will transform into the "energy representation" L 2 ([0, oo), dE) for the arrival time case, and into the "momentum 
space" L 2 (R, dk) for the other two cases. Notationally, this is conveniently implemented by defining the relevant basic 
operators (multiplication in particular) in the generic £ 2 (R), and this is done in the following. 

A. Basic notations 

For any Hilbert space Ti, we let C(H) denote the set of bounded operators on H. Let M be the multiplication 
operator acting in L 2 (R) as (Mip)(x) = xtp(x), on its domain of selfadjointness. Let D be the differential operator 
i-^:, likewise in L 2 (R). These operators will be used in different forms: in the "position representation" L 2 (R, dx), 
we will put Q :— M and P := —D; they are the standard position and momentum operators. In the "momentum 
representation" L 2 (R, dk), the multiplication operator M represents momentum, and in the "energy representation" 
L 2 ([0, oo), dE) it acts as "multiplication by energy". 

Let F £ C(L 2 (R)) be the Fourier-Plancherel operator, i.e. the unitary operator with 

{Ftp)(y) := [ e- tvx ip{x) dx, if € L 2 (R) n L X (R). 

V27T je 

The operators M and D are well-known to be connected via the operator equalities 

D = FMF* = —F*MF, (1) 

(See e.g. [l], pp. 106, 112]; it will be crucial to get the signs correctly.) We will also denote ip — Ftp, and %p = F*ip, 
for tp e L 2 (R). For any operator A in L 2 (R) (bounded or not), we denote A := FAF*, and A — F*AF. 

In the physical context, we use the Fourier operator F in the natural way as F : L 2 (R,dx) — ► L 2 (R, dk), so that 
the notations ip, ip are as usual. Also the meaning of A and A should be clear: if A is an operator in the position 
space, then A is the corresponding operator in the momentum space, and if A acts in the momentum space, then A 
is how it acts in the position space. In particular, Q = D acts as a differential operator in L 2 (R, dk), while P = M is 
the multiplication. 

For any Borel function / : R — > C, the operator f(M), as defined via the spectral calculus, is simply the multi- 
plication by /, on its domain {ip g L 2 (R) | J \f(x)tp(x)\ 2 dx < oo}. We let P+ = X[o,oo)(M) and P_ = X(-oo,o](M) 
(where \x is the characteristic function of a set X C R). Then P± £ £(L 2 (R)) are projections, and we do the 
obvious identifications P + L 2 (R) = L 2 ([0,oo)), P_L 2 (R) = L 2 {{— oo, 0]). Using the above defined notation, we have 
P± = F*P ± F. Since F 2 is the parity operator, we also have P T = I - P± = FP ± F* = P±. 

The reason for introducing these projections is that we will frequently need the subspaces of positive and negative 
momenta. In L 2 (R,dk), these are just P + L 2 (R,dk) = L 2 ([0,oo),dk) and P_L 2 (R,dk) = L 2 ((-oo,0},dk), while in 
the position space L 2 (R, dx), they are the images of the projections P±. 

B. Covariant observables 

Each of the three approaches is formulated in terms of positive operator measures, defined on the Borel c-algebra 
B(R) of the real line. We proceed to define this concept. 
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Let Ji be a Hilbert space. A set function E : B(R) — > C{Ji) is said to be a positive operator measure (POM) if E is 
strongly (or, equivalently, weakly) cr-additive, and < E(B) < I for all B G B(R). For any pair of vectors ip,ip eJi, 
and a POM E : B(R) — > we can associate the complex measure £> i— > E^ tV {B) := (ip\E(B)<p) . 

A POM will also be called observable, when a quantum system is associated with the Hilbert space Ji. The physical 
meaning of this is imported by postulating that for any state operator p (i.e. a positive operator of trace one), the 
number Ti[pE(B)] is the probability that the measurement of E yields a value from the Borel set B, given that the 
system is prepared into the state p. Note that here we do not require an observable to be normalized in the sense 
that E(R) = I. The positive operator I — E(R) is simply interpreted as corresponding to the event of "no detection". 

We will need two kinds of observables, arrival time and localization observables. In the first case the problem of 
time in quantum mechanics is obviously involved. Without delving into the long history of this question (see the 
references given in the introduction), we recall that the use of POMs is forced by the fact that there is no selfadjoint 
operator giving eligible "time" probability distributions. 

A time observable associated with the free evolution Hq := P 2 is a POM E : B(R) — > C(L 2 (R, dx)) satisfying the 
covariance condition 

e ltHo E{B)e- ltH ° = E{B + t), for all B G B(R), t G R. (2) 

This encodes the minimal requirement that the measurement of E performed at time t > gives a result from the 
range [iijtg] with the same probability as the measurement of E at t = gives a result from the shifted range 
[ti +t,t 2 + t}. 

For an arrival time observable, we additionally require that E(R) is the projection onto the subspace of positive 
momenta, i.e. E(R) = P + . This is simply because the "arrivals" are supposed to be coming only from the left, so the 
negative momentum part is not detected (see [27| for a more general formulation of screen observables.) 

In the second (and third) approach, we need localization observables. The standard localization observable is given 
by the spectral measure E® of the position operator Q. However, in order to emphasize the mathematical similarity 
of the approaches, we consider general localization observables, i.e. POMs E : B(R) — > £(L 2 (R, dx)), satisfying 
translation covariance and velocity boost invariance: 

e- up E{B)e up = E(B + t), e UQ E(B)e- UQ = E(B), for all t G R. (3) 

Such observables have been studied in the context of approximate (or imprecise) position measurements (see e.g. 
@, SBS|)- In particular, they are all known to be of the form E = p * E® , where p : B(R) — > C(Ji) is a probability 
measure, and the convolution is defined in terms of the associated complex measures. 



C. The tunnel potential and scattering in one dimension 

Having defined the covariance concepts, we move on consider the tunnel potential. Quite naturally, the essential 
quantity will turn out to be the transmission amplitude associated with the scattering from the potential. The relevant 
information will be given in Theorem Q] below. 

Let Hq := P 2 be the free Hamiltonian (in the position representation). For our purposes, a tunnel potential is a 
(measurable) function V : R — > M such that 

(i) V is compactly supported and bounded; 

(ii) The Hamiltonian H = Ho + V has no eigenvalues (e.g. V is positive). 

Condition (i) assures that the tunnel is strictly localized in some interval (xq,xi), and does not form an impenetrable 
barrier. The second condition means that it actually acts as a barrier rather than e.g. a well. In order to not to 
exclude the square barriers typically used in the context of tunneling, we have not required continuity for the potential. 

Next we need to recall some basic facts of scattering theory in one dimension. Under the above conditions defining 
the tunnel potential, the Hamiltonian H = Hq + V is selfadjoint, with purely absolutely continuous spectrum [0, oo). 
In particular, there are no bound states. The wave operators 

Q±=s- lim e ltH e-' ltH " 

t — >±oo 

exist and asymptotic completeness holds, i.e. Q + (H) — Sl-(H) — L 2 (R,dx). The operators f2± are unitary, and the 
unitary operator S = fl*Lfl, which connects incoming and outgoing asymptotics, is called the scattering operator. We 
have the following intertwining relations. 

e- ltH n± = n± e - ltHo , tel; (4) 

e ^tHo S _ Se - itHo , tER. (5) 



4 



The last equality implies that S commutes with Hq = P 2 , and this gives rise to a decomposition of S : letting P denote 
the parity operator, each of the four operators P + SP + , P + SP_P, PP_SP + , and PP-SP-P, acts on P + L 2 (R, dx), 
and commutes with the momentum P. The corresponding momentum space operators thus act multiplicatively on 
L 2 ([0, oo), dk); we will denote them by 

T r {P) := PP_SP_P, Ri{P) := PP-SP+P 

R r (P) := P+SP_P, Tt(P) := P+SP+. 

(Recall the notation: in L 2 (R, dk) the momentum operator is P — M .) The four functions Ti, T r , Ri,R r : [0, oo) — > C 
thus defined are measurable, and (essentially) bounded by one. The functions Ti and T r are called the coefficients of 
transmission, while Ri and R r are the coefficients of reflection. By denoting 

5(k):=( T ^ k>0. 
\R r (k) Tiik))' ~ 

one gets an explicit form for the action of S in the momentum space: 

fe))= sw te))' k> -°- * EiS(M) - (6) 

The fc-dependent matrix S(k), k > 0, is called the scattering matrix for H. It mixes the positive and negative 
momentum components of the "initial" asymptotically free state to produce the corresponding "final" free state, 
having "transmitted" and "reflected" parts. 

The structure of transmission and reflection coefficients is investigated via the stationary scattering theory: there 
exists, for each k <G R, two solutions <f>i(x, k) and 4>-z{x, k) of the differential equation 



d 2 

ip(x) + V(x)ip(x) = k 2 ijj(x), (7) 



dx 2 

analytically depending on k, and satisfying 

V ; \Ti(k)e lkx , x>xi ! 

[T r (k)e- ikx , x<x 
Mx,k) = ikx „ lkx . , (9) 

I R r (k)e LKX + e lKX , x > x\ 

where x$,xi G R are any two points such that the support of V is included in [xq,Xi). We refer to [IcJ for the 
stationary theory. The functions T r , Ti, and Ri, R r appearing here are exactly the transmission and reflection 
coefficients we defined via the "time- dependent" theory 2 . 

We need the following properties of the scattering matrix ficl Theorem 1]. Let C+ C C stand for the open upper 
half-plane, i.e. C+ = {lo G C | Re(w) > 0}. 

Theorem 1. Let V be a tunnel potential. 

(a) The scattering matrix S(k) is unitary for all k G K, k i— * S(k) is continuous, and we have T := T\ = T r , 
T(-k) = T(k), Rr(-k) = lUkj, Ri(-k) = R[(kj. 

(b) The transmission amplitude T can be extended to the upper half plane C+ in such a way that T is continuous 
in C+ U R, analytic in C+, and satisfies \T(k)\ < 1 for all k £ C U R. 

Remark 1. The restriction for compactly supported potentials is not really necessary for the scattering approach; 
we could just as well use a potential with no bound states, and sufficiently rapid decrease at infinity to ensure that (a) 
the Hamiltonian is a well-defined selfadjoint operator, (b) wave operators exist and are complete, (c) the stationary 
theory works (with the relations and ([9]) understood as "asymptotically" valid), and (d) the connection between 
the "time-dependent" and stationary pictures is secured. Specific conditions for each of these requirements can be 
found in standard literature (see e.g. 



pictures is 
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2 We defined T r (k), T;(fc), Ri(k), and R r (k) for positive k\ here they are extended to negative k by e.g. T;(fc) = T;(— fe). 
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For the "time zero approach", which does not directly involve scattering theory, we will use the expansion of the 
evolution e~ ltH ip in terms of the basic solutions (pi(-,k). Such an expansion is traditionally used in the context of 
stationary scattering theory; the basic solutions are called "improper eigenfunctions" of H. In general, the existence 
of the expansion is a highly nontrivial problem, which has a long history (we only mention the old work of Titchmarsh 
as well as some relatively recent papers @, @, EH)- We will need the expansion only for tunnel potentials 
(compactly supported and bounded), for which it is known to hold, according to the references just mentioned. (As 
in the case of asymptotic completeness, the problems arise mainly for slowly decaying potentials.) 

For ip belonging to the Schwartz space of rapidly decreasing functions, we define 

{k,tp)i-=— i= / (j)i{y,-k)il){v)dx = —= I (p i {y,k)ip{y)dx. (10) 
\/ 2ir Jr v 2ir Jr 

Then 

(e- itH iP)(x) = -= / -[(k^) 1 Mx,k) + (k,iP) 2 fo(x 1 k)]e- uk2 dk. (11) 
Note that the absence of bound states is reflected in this expansion. 

III. PREPARING THE "INITIAL STATE" OF THE PARTICLE BEFORE THE TUNNEL 

In the introduction we already emphasized the importance of identifying the initial states to be the same in both 
evolutions. In the first two approaches, the identification is done by means of the scattering theory; for any given 
vector state ipo € L 2 (R,dx), with ip Q € L 2 ([0, 00), dk) (i.e. positive momenta), we find ip £ L 2 (R, dx) such that 
g-itff^ ^ e ~ ltHo ip asymptotically at t — > —00, in the sense that the difference goes strongly to zero at this limit. 
This just means ip — fl-ipa. Note that here ipo is not the initial state, because the "initial time" is considered to be 
t = — 00. According to a well-known result called "scattering into cones" [12], this setup means that the particle is 
initially localized "far to the left" of the potential at t — > —00, and is "going to the right" at any time t. 

Obviously, the pure state ipo can also be replaced by a general state operator po. Then the condition of positive 
momenta is Tr[/?oP + ] — 1, or, equivalently, P+pqP+ — pq. 

In the "time zero approach", we take an interval (xq,xi) which includes the support of the tunnel potential V, and 
at the initial time t = we prepare a state ipo G L 2 ((~ 00, xq], dx). Then for t > the tunneled and freely evolved 
states are simply ipo and e~ ltH °ipo, respectively. We let Pj n it denote the projection onto L 2 ((— 00, xq], dx), so 
that we can state the initial condition for a general state po as Trfpo-Rnit] = 1- 

IV. THE "DETECTION" OF THE PARTICLE AFTER THE TUNNEL 

We describe here in detail the detection method in each of the three schemes; in each case, we end up with two 
relevant probabilities, corresponding to the tunnel particle and the free particle, respectively. 

A. Approach I: arrival time 

For an arrival time observable E, the number Tr[E((— 00, t])p] is interpreted as the probability that a particle whose 
state is p at t = arrives at a certain point (which depends on E) during the time (— 00, t]. As explained in the 
introduction, the idea is to compare the arrival time probability of a particle moving in the presence of a potential, 
with the corresponding probability of a freely evolving particle, with initial states identified as above. 

For a free particle, an arrival time observable Eq : B(M) — * C(H) must satisfy the covariance condition |(2|), and 
the additional condition that E(R) = P + . As explained before, this means that the observable is only sensitive to 
positive momenta; particles traveling "to the left" will not be detected. The correct time observable E for the evolution 
according to H should satisfy ([2]) with Hq replaced by H, because H generates the time translations for this system. 

With Eq and E chosen this way, and given a pure state ipo as in the preceding section, with ip = fl-ipo, the arrival 
probabilities to be compared are of the form (ip\E((— 00, t])ip) and (iPq\Eq((— 00, t])ipo). In order to ensure that the 
comparison is meaningful, the observables E and Eq have to be "the same after the scattering event", i.e., at large 
times t — > 00. Accordingly, we require that for any given ipo € "H, 

{(p\E(B)tp) = (<p \E (B)<p ), for all B € B(R), 
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with ip — fl+ipo. This just amounts to saying that arrival time probabilities corresponding to the two evolutions 
should coincide for states which will become asymptotically equal at t — > oo. This condition is equivalent to the 
requirement E(B) = n + E (B)£l+. Note that for any arrival time observable E corresponding to H , the observable 
B h-> E(B) := n+E (B)n* + indeed satisfies © with H replaced by H because of Q. 

Hence, in the end we actually need only the observable Eg, which satisfies we compare 

(V|S((-oo,t])V) = (Sip \E ((-oc,t])Sip ) 

with 

(V>o|£o((-oc,i])Vo}. 

Moreover, since E ((—oo,t]) < P + , and S P + L 2 (R, dx), we can simply replace S by the transmission amplitude 
P + SP + = T(P) in the former. This just means that because the observable is only sensitive to positive momenta, it 
does not "see" the reflected part of the state. Using an arbitrary state po with Tr[po-P+] = 1, we thus compare 

Pi(t) ~ Tr[T(P)*E ((-oo,t])T(P)p ], 

p?(i) :=■&[£*,((-«>, f])/*], (12) 
where the index I refers to this first approach. 



B. Approach II: localization measurement 

Here we do not have the problem of identifying the observables; we make the same localization measurement E 
for both tunnel and particle case, at a large preset time t, At this time, the corresponding states are e~ ltH ifj and 
e -itffo^ Q _ Here "large" time means that we are in the asymptotic regime, i.e. we identify e~ ttH ip ~ e~ ltH ° Stpo at 
t — ► oo, in the sense that the difference goes strongly to zero. 

Now the reflected part P-Sipo does not contribute to the localization measurement since it "moves to the left" while 
we localize in [a, oo). In the case of sharp localization (corresponding to the spectral measure of Q) this is a 
well-known consequence of the "scattering into cones" - theorem, and can be derived from the asymptotic form for 
the free propagator (see e.g. [23, p. 60]): 

s- lim E Q ([a,oo))e- ltHo P- =0. 

t — >oo 

As we mentioned when introducing the localization observables, each of them is a convolution of the sharp localization 
E® with a probability measure. Using this fact, the same limit result is easily proved also for the general case: 

Lemma 1. Let E : B(M) — ► C(L 2 (R, dx)) be an arbitrary localization observable, and a € K. Then 

s- lim E{[a l0 o))er ltHo P- = 0. 

t — ►OO 

Proof. Suppose that ip € P_L 2 (R, dx), \\p\\ = 1, and put ip t = e-' ltH «tp for t € R. Let p : B(R) -> [0, 1] be a probability 
measure, such that E = p* E® . For any ijj 6 H let E^ m denote the complex measure B t— > (tp\E(B)(pt) , and define 
E 1, Vt similarly. Then we have E^ t9t = p* E $ m , so that E^ m ([a, oo)) = (p x E^ m ){A), where A = {(x,y) € R 2 \ 
x + y > a}. Let e > 0. Since p is a finite measure, there exists an b G (0, oo) with p((— oo, —b) U (b, oo)) < |. Let 
to > be such that for t > to we have ||i? Q ([a — b, oo))ipt\\ < § • Then for t > to, and any unit vector ip £ H, we have 
|£ , ^ v J([a — b, oo)) < | (where the | • | stands for the total variation of the measure), and consequently, 

\E^ Vt ([a,oo))\ < (px\Elj)(A) 

< / /((-oo,-&)U(6 ) oo))|^ V3 J(R)+M(R)|^ t |([a-&,oo))<e. 
This implies ||J5([a, oo))ip t \\ < e, so the proof is complete. □ 
Thus, we may identify 

E([a, oo))e- itH tp - E([a,oo))e- ltHo S^J Q - E([a,oo))e^ tH "T(Pyjo, 

in the sense of strong asymptotic convergence. Hence, in this second approach, we want to compare the localization 
probabilities 

p u (a) := Tr[T{PyE{[a,^))T{P)e-* tH ° P() e HH °]; 

p° u (a) := T r {E({a,^))e- ttH °p Q e ttH °], (13) 
for arbitrary states po with Tr[/?o-F+] = 1. 
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C. Approach III: time zero 



Here we make a sharp localization measurement, with the observable , at time t > corresponding to an interval 
[0,00), where a > X\. Recall that t = is the initial time, the support of the potential is contained in {xq,X\), and 
the initial state po satisfies Tr[poPnit] = 1, where P n i t = E^((~ 00, xq]). In view of the result ((15)1 below, it will be 
convenient to put P nna i := E^([xi,oo)), 

If po = \i/jo)('4>o\, where tpo belongs to the Schwartz space, we can use the expansion (fTT]l . Since the support of i/jq 
is contained in (—00, xq], we can express (V>o, k)i in terms of the Fourier transform; indeed, by and ([9]), as well as 
Theorem Q] (a) , we get 



(ipo,k)i = —}= I ip (x)(j)i{x,~k)dx = tpo(k) + Ri(-k)tp (-k), 

' Z7T 



/Xq 
ipo(x)(j)2{x,-k)dx = T(-k)ip (-k). 
-00 

The expansion (fTTjl now takes the form 

(e- ltH ij )(x) = -= / -{(Mk) +Ri(-k)M~k))Mx,k) + T(~k)M-k)M^, k)} e - uk2 dk. 
Vztt J -00 -a 

Using again the relations iJH) and ® we get, for x>X\, 

{e~ UH ^){x) = -= / -[^ (k)T(k)e* kx + M~k)T(~k)e-* kx 
+ [Ri(-k)T(k)+T(-k)R r (k)]tPo(-k)e lkx ]e- ltk2 dk 

1 r°° i 

' ~[Mk)T{ky kx + $ {-k)T{-k)e- u > 



2tt i_oc 2 



where the second equality is obtained by Theorem [TJ (a) . Hence, in this case, the evolution is simply given by 

1 r°° 

(e- UH ip )(x) = -= / e lkx T(k)-ip {k)e- uk2 dk, x > x u (14) 

V27T J -00 

or, equivalently, 

Pfinaie- Itff ^ = Pan,iF*T{P)e- up2 FiP . 

Since the operators F* , T(P), e~ ltp , and F are all bounded, and the Schwartz space functions with support contained 
in (— 00, ato] are dense in L 2 ((—oo,xq]), the above relation is equivalent to the operator equality Pa n ede~' ltH Pnit = 
P tLn3 xF*T{P)e- itM2 FP init , i.e. 

Pfinaie-^Pinit = PfinalTXPK^Pinit. (15) 

We want to compare the sharp localization probabilities pm(a) := Tr[E^([a, oo))e~' ltH poe ltH ] and pj n (a) := 
Tr[P Q ([a, oo))e" itifo poe l * ff °], where Tr[p Pmit] = 1, and t > 0, a > x x are fixed. Using ([Tg|. and noting that 
finitPo-Pinit = Po, PfinaiP ([a, oo)) Pfinai = P Q ([a, oo)), we immediately get 

pin(a) = Tr[T(P)*P«([a,oo))T(P) e -" ff Voe ,;tffo ]; 

p° n (a) = Tr[P«([a, ocOJe-^poe**"]. (16) 

Note that although this expression is similar to the one in the above scattering theory with localization - approach, 
here po has also negative momentum components, and thus also the values of T for negative argument are used. 



V. A MATHEMATICAL THEOREM 



By looking at the probabilities lfl"2|). (fl3|l . and lfl6|) it is clear that the problem of comparison is similar in each case. 
As we have already pointed out, we only need to transform the covariant observable in each case into the spectral 
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representation of the associated generator, which is Hq = P 2 in the arrival time case, and P in the localization case. 
Since the spectrum of H is [0, oo), while the spectrum of P is (—00,00), we need to consider both L 2 ([0, 00)) and 
L 2 (R). In their respective spectral representations, the operators H and P act as the multiplication operator M. 

In order to deal with both cases, we formulate the essential result (Theorem [2] below) for a subset I G B(M) (either 
R or [0,oo), in practise), and the corresponding subspace L 2 (I) = xi(M)L 2 (M). For a Borel function / : R — ► C, this 
subspace is obviously invariant under the multiplication operator f(M) for any / : R — > C, this operator being the 
multiplication by the restriction of / to I. For notational simplicity, we will use the symbol f(M) also to denote this 
restricted operator. Such a restriction will be needed for the transmission amplitude, when we restrict it to positive 
momenta. 

It should be clear from the above discussion that here we do not fix the physical interpretation of M; consequently, 
also the operators P± will appear in the following lemma without such interpretation. Since these projections are not 
involved in the statement of Theorem [2] (which is the only thing we need to refer to) they can just be regarded as 
mathematical auxiliaries in this section. 

The following lemma contains the essential mathematical ingredient we need. The symbol C+ C C stands for the 
upper open half-plane., i.e. C+ = {uj G C | Re(cj) > 0}. 

Lemma 2. Let g : R — > C be a measurable function, and suppose that g : C+ U R — > C is an extension of g such that 
(i) g is analytic on C+; 
(ii) \g(uj)\ < 1 for all lo G C+ U R; 
(Hi) lim;,^o+ g{ a + ib) — g{a) for almost all a G R. 
Then 

g{M)*P + g{M) < P+. 

Proof. Put Ji = L 2 (R) for convenience. We first note that it is sufficient to prove 

g(M)PJH C P-H, (17) 

Indeed, QU implies (ip\g{M)* P+g(M)ip) = (g(M)ip\P+g(M)<tp) = for tp G P+iU) 1 - = P-(H), ip G H, that is, 
g(M)*P+g(M)H C P+H; this gives g(M)*P + g(M) < P + , because g(M)*P+g(M) is a positive operator with norm 
less than one by (ii). 

Now ip G P-(H) if and only if F*ip is supported in (—00, 0], or equivalently, if and only if Ftp is supported in [0, 00). 
To prove (fl7f . let 93 G P-(H), and define (p : C+ — > C via 

-L f e^(Fip)(y)dy, lo G C+, 

where the integral exists because y 1— > \e lujy \ = e~ Im (^ v is now square integrable over [0, 00). Then (p G H + , where 

H + = < h : C+ — *■ C I h analytic, sup / \h(a + ib)\ 2 da < 00 > 
L 6>o Jr ) 

is called the Hardy class; see e.g. fljjl p. 161-162]. According to this reference, the elements h G H + are characterized 
as precisely those functions h : C+ — * C for which there exists a function / G L 2 (M.) supported on [0,oo), such that 

Here / is recovered via the L 2 (R)-limit F*f = limh^o hb, where hb(a) := h(a + ib). Since (p G H + , the assumptions 
(i) and (ii) clearly imply that also lj ^ g(u)(p{uj) is an element of H + . Hence, there exists a tp G P—CH) with 

§(w)fp(w) = -i / e^iF^dy, u>eC+. 

To conclude the proof of (fl7| . we have to show that V = g(M)<p. According to the definition of H + , we have 
Tpbi^b G L 2 (R) for any 6 > 0, where ipb( a ) '■= g( a + ib)<p(a + ib) and <^fc(a) := tp(a + ib). Let g;, be the function 
a h- > g(a + i6) for each 6 > 0, so that V'fc = 9b(M)(pb, for 6 > 0. As mentioned above, we then have the L 2 (R)-limits 
lim^o ipb = ip and lim b ^ fb = V- Since the family of multiplication operators {gb(M) | b > 0} is uniformly bounded 
because of (ii), and gb(M) tends strongly to g(M) as b — > (by (iii), (ii) and the dominated convergence), it follows 
that ip = g{M)tp. The proof is complete. □ 
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Next we need information on the structure of covariant observables. This is well-known (see e.g. [171. l27l|). For the 
purposes of this paper, it is convenient to state it in the following form, which can be immediately specialized from 
the general construction procedure, obtained by combining Mackey's imprimitivity theorem with a dilation argument 
& 

Lemma 3. Let I C R be a Borel set, and make the usual identification L 2 (I) — %/(M)L 2 (R) C L 2 (R). Let 
E : B(M) — > C(L 2 (I)) be a positive operator measure with the covariance property 

e~ laM E{B)e iaM = E(B + a), for all B e B(R), a e R. 

Then E is of the form 

E{B) = V(Fxb(M)F*), for all B € B(WL), 
where W : £(L 2 (R)) — > C(L 2 (I)) is a positive linear map satisfying 

*(/(M)*A/(M)) = f(M)*V(A)f(M), 

for any bounded Borel function f : R — * C and A 6 £(L 2 (R)). 

With this result, the ingredient given by Lemma [2] can be imported into the covariant setting: 

Theorem 2. Let I € B(R), and let E : B(R) — > C(L 2 (I)) be a POM with the covariance property 

e~ laM E(B)e laM = E(B + a), Be B(R), a e R. 

Lei <? : R — > C satisfy the conditions of Lemma® Then 

g(M)*E([a, oo))g(M) < E([a, oo)), a e R. 

Proof. First note that by covariance, £/([a,oo)) is unitarily equivalent with £7([0, oo)) via the unitary e wM , which 
commutes with g(M). Hence, it suffices to prove the inequality for t = 0. To this end, let \& : £(L 2 (R))) — > L(L 2 (I)) 
be the map given by the above lemma, corresponding to E. Since g is a bounded Borel function, and * is positive 
and linear, Lemma [2] gives 

<7(M)*£([0,oo))<7(M) = g(M)**(FP + F*)g(M) = *( 5 (M)*P +5 (M)) 
< =£([0,oo)), 

and the proof is finished. □ 

Remark 2. Notice that although the range of all the E(B) are here asserted to be in L 2 (I), and, consequently, only 
the restriction g\i appears in the inequality of the above theorem, it is still necessary to have g as a function on the 
whole R. This is because otherwise we could not move g{M) inside the argument of ^, and apply Lemma [21 



VI. RESULTS 



We can now apply Theorem [2] of the preceding section to the three relevant cases, in order to compare the proba- 
bilities in lfT2j) . fl3l ). and [fl6j) . 



A. Approach I: arrival time 

We pass from "position representation" to the "energy representation" by means of the unitary operator U+F, where 
the unitary U + : L 2 ([0, oo), dk) — > L 2 ([0, oo), dE) is given by 

(U + ip)(E) = 2-?E-itp(y/E), for all tp € L 2 ([0, oo), dk). 

The transmission amplitude T(P) — T(M) (which already acts in L 2 ([0, oo), dk)), is transformed into U+T(M)U+ = 
T'(M), where V : R -> C is given as T'(E) = T(y/E), and the arrival time observable E : B(R) £{L 2 (R,dx)) 
transforms into the observable E' : B(R) — * C(L 2 ([0, oo), dE)), where 

E' Q (B) = U+FE (B)F*U* 
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Note that since the range of Eq(B) is included in P + L 2 (R, dx) (positive momenta), the range of FEq(B) is included 
in £- 2 ([0, oo), dk), so that this is indeed well-defined. The probabilities in ffl2|) now have the form 

pi(t) = Tr[T'(M)*E' ((-oo,t])T'(M)U + p Ul], 
p°(t) = Tr[E' ((-oc,t])U + p Ut}- 

Choosing the square root branch 

{re ie | r > 0, -vr < 9 < tt} 3 re w h-> ^Fe^ e € {re l9 \ r > 0, -tt/2 < < vr/2}, 

and using Theorem [TJ (b), we see that T" = T(y/^) satisfies the conditions of Lemma [21 Due to the intertwining 
relations U + e ltM — e ltM U + , Fe ltH ° — e ltM F (which hold for all t S M), the observable E' Q satisfies the covariance 
condition 

e ltM E' (B)e~ ltM = E' a {B + t), t e K, 

which implies that B i— > E' (—B) (and not .Eg itself) satisfies the covariance condition of Theorem [51 This reflection 
simply means that the conclusion of Theorem [2j holds for observable E' Q with the intervals — [— t, oo) = (— oo, t], and 
we get 

pi(t) < p?(t), i e R, (18) 

for any state po with Tr[po-P+] = 1- This means that the probability of arrival by the time t is never larger for the 
tunneled particle. 

B. Approach II: localization measurement 

Here the application of Theorem [JJ is more straightforward: In comparing the probabilities lfT3j) . we only need 
to pass to the momentum representation; define E : B(M) — > C(L 2 (R, dk)) as E{B) = FE(B)F*. Because of the 
intertwining Fe~ ltp = e~ ltM F, this observable now satisfies the covariance condition of Theorem [JJ In addition, T 
satisfies the conditions of Lemma [2] by Theorem Q] (b) , so Theorem [2j immediately gives 

T(M)*E([a,oo))T(M) < E([a,oo)), for all a e M. (19) 
Applied to the probabilities l|13p . this implies 

Pn(a) < p?i(a), a G R, 

for any state po with Tr[po-P+] = 1- (Recall also that the localization measurement is made at large time t; for small 
t, the comparison is just between the transmitted parts.) The result means that the probability of having passed the 
point a at time t (in the sense of the particular localization used) is never larger for the tunneled particle. 

C. Approach III: time zero 

Finally, consider the probabilities ifTB]) . Since E® is a localization observable, the inequality (fT9| immediately 
applies also to this case, and we get 

Pm(a) < Pin(a), a>xi, 

for any state po with Trfpo-Pinit] = 1- This has the same meaning as in the above case, except that here the localization 
is only understood in the sharp sense, and the time t at which the measurement is performed can be any t > 0. 
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